Newton's cradle and entanglement transport in a flexible Rydberg chain 
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In a regular, flexible chain of Rydberg atoms, a single electronic excitation localizes on two atoms 
that are in closer mutual proximity than all others. We show how the interplay between excitonic 
and atomic motion causes electronic excitation and diatomic proximity to propagate through the 
Rydberg chain as a combined pulse. In this manner entanglement is transferred adiabatically along 
the chain, reminiscent of momentum transfer in Newton's cradle. 

PACS numbers: 32.80.Ee, 82.20.Rp, 34.20.Cf 



Nearly loss-loss transfer of momentum and energy 
through a linear ehain of masses was studied as early as 
in the 17"^ century, exemplified by Newton's cradle 
Newton's cradle, the well known table-top experiment of 
popular physics, in which the first and last member of a 
suspended chain of metal balls periodically exchange os- 
cillation energy, has recently found a modern implemen- 
tation with ultra-cold atoms ■ Stimulated by research 
on photosynthesis Q and organic dye aggregates an- 
other type of energy transfer has received more recent in- 
terest, namely the propagation of internal excitation due 
to electromagnetic interactions. Experiments indicate ro- 
bust, coherent energy transport of this kind in systems 
as large as photosynthetic light harvesting complexes Q , 
opening up the avenue to consider the propagation of a 
particularly "fragile" quantum phenomenon: entangle- 
ment 0] ■ Here we show for the example of a chain of 
ultra-cold Rydberg atoms that mechanical momentum 
transfer interlinked with coherent excitation migration 
can result in efficient transport of entanglement. 

Rydberg atoms have recently received much attention, 
to a large part due to their strong long-range interactions, 
with diverse co nseq uences from dipole-blockade [8| and 
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anti-blockade 0, llOj over long range molecules [UJ, 
classical motion due to Van-der-Waals interactions \lc 
In contrast to the latter, resonant dipole-dipole interac- 
tions [H - 17 1 intimately link motion and excitation trans- 
port. Within an essential state picture, where only two 
electronic Rydberg states per atom, labelled | a ) and \b), 
are taken into account, the transfer of excitation can 
be adequately described by using the exciton theory of 
Frenkel Q . For a pair of atoms separated by a distance R, 
dipole-dipole interactions have a Hamiltonian with struc- 
ture H = V{R){\ab){ba\ + \ba){ab\), where V{R) scales 
like R~^. It describes electronic excitation transfer, since 
a transition of the first atom from | a ) to \b) is accompa- 
nied by the reverse transition of the second atom. To see 
how this also induces mechanical forces, we consider a su- 
perposition of two-atom states like \ab)±\ba), which are 
excitonic eigenstates of H with eigenvalues ztV{R) that 
parametrically depend on the distance R. These provide 
adiabatic potentials for the nuclear motion. The charac- 



ter of the motion (eg. fully repulsive or fully attractive) 
depends on the exciton state 17[- Adiabatic motion of 
atoms in a longer chain preserves the exciton character. 
Since the exciton state for more than two atoms depends 
on the atomic positions, excitation transport and motion 
become interlinked. 

In detail, we study the effect of resonant dipole-dipole 
interactions on a regular linear chain of Rydberg atoms. 
Initially we impose a perturbation in the distances be- 
tween the atoms by placing two atoms close together, 
with a localized exciton state built on this diatomic prox- 
imity, choosen repulsively. We demonstrate a strong cor- 
relation between the resulting exciton dynamics and the 
motion of the atoms: The combined pulse of atomic dis- 
placements and localized electronic excitation propagates 
adiabatically through the chain in a manner reminiscent 
of Newton's cradle. 

We treat this complex many-body problem using a 
mixed quantum-classical approach (Tullys surface hop- 
ping method 3, 31). It allows us to determine the dy- 
namics of the atomic wave-packets together with the elec- 
tronic excitation transport in order to quantify the en- 
tanglement in time. For short chains, where a full quan- 
tum mechanical treatment is possible, the numerically 
exact solution is in perfect agreement with the quantum- 
classical result. 

We study a linear chain of N identical atoms with mass 
M and denote by i?„ the position of the nth atom (nu- 
clear coordinates). All but one of the N atoms shall 
be in a Rydberg state |i^s), i.e. with principal quantum 
number ly and angular momentum 1 = 0. Just a sin- 
gle atom is in the state {vp), i.e. has angular momen- 
tum 1 = 1. The latter will be called the "excited" state 
hereafter. It can migrate along the chain by means of 
dipole-dipole interactions [3], which conserve the num- 
ber of excitations. We can restrict ourselves to the single- 
excitation Hilbert space, whose electronic part is spanned 
by 1 7r„ ) = \ s- ■ -p- ■ ■ s), see Fig. [T] (b), since for the sce- 
nario we consider, transitions to other states are negligi- 
ble. The distance Rnm = \Rm — Rn\ between the atoms 
n and m is so large that the overlap of their electronic 
wave functions can be neglected. The total Hamiltonian 
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of the system is 



TV 



(1) 



where R = 
tions. The electronic Hamihonian 



. , Rjy)'^ is the vector of nuclear posi- 



(2) 



contains the dipole-dipole coupling between atoms n 
and m. We consider the case with all atoms in 
mi = azimuthal angular momentum states, such that 
Vnm{Rnm) = — P-^ / Rnm without angular dependence [l6| . 

Our numerical calculations use an atomic mass M = 
11000 au (which is roughly the mass of Lithium) and a 
transition dipole moment /i = 1000 au, corresponding to 
transitions between s and p states with n sa 30 . . . 40. 

The full many-body problem posed by the Hamilto- 
nian ([T]) becomes quickly intractable as the number of 
atoms N is increased. For small N however, it is no prob- 
lem to directly solve the equation of motion. Expanding 
the full wave function in electronic (diabatic) states ac- 
cording to I ^'(R) ) = 4'ni^)\ TTn ), wc arrive at the 
Schrodinger equation (in atomic units) 
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In order to validate the semi-classical method pre- 
sented below, which in turn will be faithfully used 
for longer chains, we solve Eq. ([3]) for iV = 3. In 
our figures we will not show the full iV-dimensional 
nuclear wave function but focus on the more intu- 
itive total atomic density, which is given by n{R) = 
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denotes integration over all but the j — th nuclear co- 
ordinate. The density n{R) gives the probability to find 
an atom at position R. 

The diabatic representation of the wave func- 
tion allows a straight forward propagation for short 
chains. For longer chains and for the interpretation 
of the results, the adiabatic representation |5'(R)) = 
J2n=i 't'n{'R) \ <y5n(R) ) IS hclpful. Here the adiabatic basis 
|^„) is defined via 7?^i(R)| (^„(R) ) ^ [/„(R)| y5„(R) ). 
For each R there are N excitonic eigenstates | ipn (R) ) 
labeled by the index n. The corresponding eigenen- 
ergies f/n(R) define the adiabatic potential surfaces. 
The two representations are related by (/)„(R) = 
E,„(¥'n(R)|7r,„)0,„(R). 

For long chains, we solve the time-dependent 
Schrodinger equation with Hamiltonian ([T]) using a mixed 
quant um/ classical method, TuUy's surface hopping algo- 
rithm 18, 1^. In this approach an ensemble of trajecto- 
ries is propagated, and each trajectory moves classically 
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FIG. 1: (color- online) (a) Sketch of the initial total density 
distribution of A'^ = 5 Rydberg atoms, (b) Visualization of the 
electronic state | tti ). (c) Trapping of the electronic excitation 
in the repulsive exciton state by a perturbation of the regular 
chain. Shown are the populations p„ = |(7r„ | </3rcp(R))|^ as 
a function of a/xo- 



on a single adiabatic surface J7m(R), except for the pos- 
sibility of instantaneous switches among the adiabatic 
states. 

The equations of motion read 
d ^ 

i-^Ck = Uk(JV)Ck-iY'R-dkqCq, (4) 

q=l 

Am = - Vr {ipra (R) I (R) I (R)) • (5) 



The N complex amplitudes Ck define the electronic 
state via |^'(R, i)) = J2n=i^n{t)\'fn{'R,t)) and the 
dfcq are non-adiabatic coupling vectors. Besides their 
appearance in Eq. (|4]), they also control the likelihood 
of stochastic jumps from the current surface m to an- 
other surface m' in Eq. ([5]), which is proportional to 
|R • dm'm\- Further details about this scheme can be 
found in Refs. [13, [3, [l^l- The underlying quantum- 



classical correspondence is discussed in Refs. [211,123. 
randomize the initial classical positions and velocities for 
the trajectories according to the Wigner distribution of 
the initial state, described in Eq. ([5]). This is essential 
for a correct description. 

Initially we assume that the Rydberg atoms are ar- 
ranged in a straight line. The distance between the first 
two atoms is denoted by a and assumed shorter than the 
equal distances (xq) between the other atoms, as sketched 
in Fig. m 

For fixed classical positions, one of the N eigenstates of 
the electronic Hamiltonian ([2|) leads to a situation where 



initially all atoms repel each other [17[ . In the following 
we will focus on this state, which we label with "rep". 
Since the dipole-dipole interaction between the first two 
atoms is much stronger than between all others, the exci- 
tation in this repulsive state is mainly localized on these 
two atoms. For a xq this initial state can be approx- 
imately written as (| tti ) — 1 7r2 ))/y/2. In Fig. [T] (c) the 
electronic population on the various atoms is shown as 
a function of a/xo- For our simulations, we have taken 
Xq = 5/im and a = 2/im, i.e. a/xo = 0.4. Then the 
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FIG. 2: (color online) Nuclear dynamics in the case A'^ = 3. 
The time evolution of the total atomic density n{x, t) (a) is 
shown together with a comparison of TuUy's surface hopping 
calculations (black solid) with the full quantum evolution (red 
dashed) in the other panels, (b) Spatial slice n{x, to), with to 
as indicated by the first vertical white lines in (a). Arbitrary 
units, (c) Relative population n2 = / dR|02(R)|^ (n2 = |c2|^ 
in the TuUy algorithm) on the adiabatic surface (index 2) that 
is energetically nearest to the initial repulsive one. This is a 
measure of the propensity of non-adiabatic transitions. The 
deviation of curves in (b,c) is shown magnified in the insets, 
(d) Initial repulsive state. 

dipole-dipole interaction between the last two atoms is 
about 5 times larger than for the rest of the chain and 
more than 90% of the excitation is localised on the first 
two atoms. Prom the point-like (classical) nuclear posi- 
tions, we now move to a quantum nuclear wave function. 

The spatial wave function of each atom is assumed 
Gaussian with a standard deviation cto- Hence, we take 
the complete initial wave function (i.e. containing nuclear 
and excitonic degrees of freedom) as 

N 
n=\ 

0G(i?„) - AAcxp(-[i?„ - i?o„]V2CTo), (6) 

where i?ori is the center of mass of the n-th Gaussian and 
M a normalization factor. 

To confirm the applicability of the quantum-classical 
numerical treatment, we consider the smallest non-trivial 
chain N = 3. Figure [2] shows the quantum mechanical 
probability to find an atom at a certain position, in per- 
fect agreement with the corresponding graph obtained 
with the quantum-classical hybrid approach. 

The excellent agreement between the two disparate 
methods for = 3 gives confidence that TuUy's surface 
hopping produces reliable results also for longer chains, 
such as A'^ = 7, which we consider now. The correspond- 
ing atomic motion and excitation transfer, when start- 
ing in the exciton state with highest energy are shown 
in Fig. [31 Let us first consider the atomic motion. As 
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FIG. 3: (color online) Dynamics of atomic motion and exci- 
tation transfer, (a) Total atomic density averaged over IC 
realizations. We actually plot ^n. (b) Mean trajectories of 
the individual atoms (white) and electronic excitation proba- 
bilities (diabatic populations) |cmp, Cm = J^m^nrnCm- The 
latter are encoded as the width of the copper shading sur- 
rounding each trajectory, (c) Population on the adiabatic 
surface rep (black) and individual diabatic populations, (d) 
Purity P —Tr[a^] of the reduced electronic density matrix 
a (green) and bipartite entanglement i5„,„+i for neighboring 
atoms as defined in the text. The dotted line indicates 1. 



expected, initially the two excited atoms strongly repel 
each other. When atom 2 has approached atom 3, the 
main repulsion is now between those two, causing atom 
2 to slow down and atom 3 to accelerate. In this way 
the initial momentum is transferred through the chain 
to atom 7; realizing a microscopic version of Newton's 
cradle. 

To lay the basis for the treatment of entanglement dy- 
namics, we next discuss the excitation transfer, shown 
in Fig. [3] (c), which is strongly coupled to the atomic 
motion as can be seen in Fig. [3] (b) . The excitation gets 
transferred remaining always localized on the two instan- 
taneously nearest atoms, in accordance with the struc- 
ture of exciton eigenstates outlined in [13. After 5.5 ^s 
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the momentum transferred through the chain kicks out 
the last atom, and a well defined close proximity pair 
no longer exists. The exciton state then assumes the 
shape for an equidistant chain, de-localized over the en- 
tire chain (consisting of the remaining atoms), which 
subsequently slowly spreads out. From the occupation of 
the initially populated repulsive adiabatic state, which is 
also shown in Fig.[3](b), we can deduce that the evolution 
is largely adiabatic. Non-adiabatic effects occur only in 
a small interval around t = 5.5/^s, when the last atom 
leaves the chain. 

So far we have described transfer of momentum and 
kinetic energy through the Rydberg chain, both of which 
would also occur in a classical Newton's cradle. How- 
ever, the microscopic excitation migration leads also to a 
transport of entanglement, which has no classical equiv- 
alent. In the spirit of Newton's cradle we focus on 
electronic entanglement between two atoms, which we 
will quan tify with the bipartite entanglement of forma- 



tion 
sity 



23, 24 



To this end we need the reduced den- 
= J2n.m ^rim 1 7r„ ) ( TTm | , describing the 



matrix a = Z^„,„i o-nm| 7r„ 
electronic state of the system, after tracing over the 
atomic positions. The matrix elements are given 



j(R)(/),„(R) for the full quantum 
for the surface hopping 



In the latter case denotes the trajectory av- 

T ~ 



by Cr„m = / C'^R- 

calculations and a. 
method 

erage and c„ = Cnm^m are the coefficients in the 
diabatic basis. From a we then construct the binary 
reduced electronic density matrix of atoms a and b by 
= Tr^'^^''^ [a] . The symbol Tr^'^'''^ [ • • • ] denotes the 
trace over the electronic states for all atoms other than 
a, b. The remaining reduced subspacc of atoms a and 
b is spanned by \pp), | ) , \ sp)i I ) . Due to the 
structure of the 1 7r„ ) , the only non- vanishing matrix ele- 
ments of Pab are {pS \Pab\ps) = (Taa, {sp\^ab\ Sp) = (Tbb, 
{ps\Pab\sp) = {sp\j3ab\ps)* = (7ba and {ss\Pab\ss) = 



^c^{a b} '^cc- From the matrix j3 we can derive the con- 
currence Cab = 2|crab| as Outlined in [23| . The concur- 
rence is already a measure of entanglement with < 
Cab < 1, from which we finally obtain the entanglement 
of formation < Eab{Cab) < 1 as described in |23l. |24|. 

As can be seen in Fig. [3] (d), the initially perfect en- 
tanglement between atom 1 and 2 is transported through 
the chain with only minor losses up to the point where 
the final atom leaves the chain around t = 5.5^s [isl ]. 
At that moment, the exciton state de-localizes over the 
entire chain. Entanglement is then shared among all re- 
maining atoms, with a resulting drop of bipartite entan- 
glement. 

In summary, for an aggregate of Rydberg atoms in 
form of a linear chain, we have identified a dynamical 
mode that links the motion of the atoms intimately with 
the coherent propagation of a single electronic excitation 
along the chain. Adiabatic transport ensures that the ex- 
citation remains spatially localized near a diatomic prox- 
imity passing through the chain. Akin to the transfer of 



the almost macroscopic quantities, energy and momen- 
tum, in Newton's cradle, the mode transports localized 
coherent excitation and electronic entanglement along 
the chain. For the underlying quantum many-body prob- 
lem we have demonstrated the applicability of TuUy's 
surface hopping method 3; 3 by comparison with ex- 
act calculations. 
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